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We extend the recursion formula for matrix Bessel functions, which we obtained previously, 
to superspace. It is sufficient to do this for the unitary orthosymplectic supergroup. By direct 
computations, we show that fairly explicit results can be obtained, at least up to dimension 8x8 
for the supermatrices. Since we introduce a new technique, we discuss various of its aspects in some 
detail. 
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I. INTRODUCTION 

o ■ 

In a previous work, we studied properties of matrix Bessel functions in ordinary spaceO. Here, we generalize these 
^vq ' investigations to superspace. For the introductory remarks and the mathematical and physical background relevant 
for the ordinary space, and also relevant as the basis for tha,present study, we refer the reader to Ref£3. . .. . 

In mathematics, supersymmetry was pioneered by Berezincl and, in particular group theoretical aspects, by Kacr 9 lr°l. 
The theorjj-]Of non-linear a models in spaces of superm atcix fields was developed in physics of disordered systems 
by Efetovtffl. Vepbaarschot, Weidenmiiller and ZirnbauercSnH used his approach to study mndels in Random Matrix 
Theory. In RefJj, the first .supersymmetric generalization of the Itzykson-Zuber integralEEl was given. In Ref.Ej, 
Gelfand-Tzetlin coordinates!! were constructed for the unitary supergroup. Extending Shatashvili'st3 method, the 
supersymmetric Itzykson-Zuber integral was also rederived in Ref.113 in its most general form. Using the techniques 
of RefB, such a calculation was also performed in Ref.El. 
^ — , ■ From a mathematical viewpoint, Efetov's worku is the basis for a harmonic analysis in certain supersymmetric 
coset spaces, the Efetov spaces, which are relevant for the non-linear a models. In the full superspaces, a technique 
involving convolution integrals and ingredients of the corresponding harmonic analysis was introduced in RefJlj. In 
the Efetov spaces, the theory of harmonic analysis, in both its,mathematical and physical aspects, was developed 
by ZirnbauerEa and was applied to disordered systems in Refs.e2rE3. In the present contribution, we do not focus on 
the Efetov spaces, rather we address the full supergroup spaces. The supersymmetric Itzykson-Zuber integrals and 
its application in Ref.E.3 is the simplest example of a supermatrix Bessel function appearing in this kind of harmonic 
analysis. 

The matrix Bessel functions in superspace find direct application in Random Matrix Theory. For general reviews, 
see Refs.E3EjO. In Ref.Eil it was shown that they are the kernels for the supersymmetric analogue of Dyson's Brownian 
Motion. 

The paper is organized as follows: In Sec. |l[ we introduce the supermatxix Bessel functions and collect basic 
definitions and notations. In Sec. [II, we extend the recursion formula of Ref.E3 to superspaces. Since it is one of our 



goals to demonstrate that explicit results for supermatrix Bessel functions can indeed be obtained, we present, in some 
detail, such calculations for certain supermatrix Bessel functions in Sees. |y| and |v|, respect ively . The asymptotics 



and the normalizatio n arc discussed in Sec. VI. We briefly comment on applications in Sec. VII and we summarize 



and conclude in Sec. |VIII| . Various calculations are shifted to the appendix. 



II. SUPERMATRIX BESSEL FUNCTIONS 

Similar to ordinary spacesEl, the superunitary aase. i.e. integration over the supergroup U(ki/k2), is the simplest 
one. As it was already discussed in detail in Refs.ErtJ, we refrain from reconsidering it here. Thus, it turns out that 
we may restrict Q ucse lves to the supermatrix Bessel function of the unitary orthosymplectic group UOSp{ki/2k2). As 
discussed by Ka<xac3, the supergroups U {k\/k2) and UOSp{k\/2k2) exhaust almost all classical compact supergroups, 
apart from some exotic exceptions which are of little relevance for applications. Thus, the integral we have to deal 
with is given by 
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$fei2fe 2 (s,r) 



u£UOSp(k 1 /2k 2 ) 



exp(itrgu 1 sur^d/i(u) , 



(2.1) 



distinguish the two possible forms. We will return to this issue. The matrices si, s 2 and r%, r 2 are given by 



where dfx(u) is the invariant measure. The arguments of the function (2J1 are the diagonal matrices s — 
diag (y/csi, \J — cs 2 ) and r = diag (y/cri, \J-cr-i). Here, we use Wegner's notationEj and introduce the label c = ±1 to 



(2.2) 



s 1 = diag(sn,s 2 i, . 
r x = diag(rn,r 2 i, . 



iSfcii)i s 2 = diag(si 2 l 2 , • 
,r kl i), r 2 = diag (r i2 l 2 , . 



! Sfc 22 l 2 ) , 

,r k2 2h) ■ 



There is a twofold degeneracy in s 2 and r 2 , because the matrix u 1 su or, equivalently, uru 1 has to be a real Hermitean 
supermatrixES of the form 



r(A) 



—OjiHSd) 



±1 



(2.3) 



The matrices and a^ HSd ^ have ordinary commuting entries, i.e. bosons, they are real symmetric and Hermitean 
self-dual, respectively. The matrix has anticommuting or Grassmann entries, i.e. fermions, and is of the form 



■ j CT fci J 



7 {A) 
J li 

(A)* 



k 2 i 
(A)* 
T k 2 i 



(2.4) 



We can now appreciate the meaning of the parameter c which enters the definition ( J2 .3| ) of the real Hermitean matrices. 
For c — 1, it yields the real symmetric and for c = — 1 the Hermitean self-dual matrix as boson-boson block, and 
vice versa for the fermion-fermion block. In the framework of Random Matrix Theory, we find the supermatrices 
corresponding to the Gaussian Orthogonal Ensemble (GOE) for c = +1 and those for the Gaussian Symplectic 
Ensemble (GSE) for c = -1. 

The infinitesimal volume element is given by 



fel k 2 



ki 



k 2 



m = n n n < 5 n ^ n ^ Sd) \ n d ^ sd) 
i=ij=i 



(2.5) 



The supermatrix Bessel functions ( |2.lD are eigenfunctions of a wave equation in the curved space of the eigenvalues 
s or r. As in the ordinary case, a supermatrix gradient d/da is introduced and the Laplace operator is defined by 



where d[a^ Sd " > ] is the product of the differentials of all independent elements of the quaternion a VJ 



(HSd) 



trg 



d_ 



The plane waves exp(itrgcrp) are the eigenfunctions, i.e. we have 

Aexp(itrgcrp) = — trg/9 2 exp(itrgcrp) 



(2.6) 



(2.7) 



Here, both a and p are real Hermitean. As in ordinary space, the supermatrix Bessel functions are obtained by 
averaging over the angular coordinates, i.e. over the diagonalizing group . The Laplacean commutes with the average 
and we arrive at the differential equation 



A s $ fel2 fc 2 (s,r) = -trgr 2 $ fel2 fc 2 (s,r) , 
where the radial part of of the Laplacean ( |2.6|) reads 

A., = 




(2. 



(2.9) 
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The Jacobian or Berezinian is given b; 

(1) _ \A kl ( Sl )\Ai 2 ( 1S2 ) _ jAfaCM^) 

n P =i n g =i( s Pi - is ?2) 2 n P =i n,=i(*s P i - s q2 y 

One easily convinces oneself that A s depends on c only through a factor y/c. Thus, without loss of generality, we set 
c = 1 and omit the index c. 

At this point, an important comment is in order. The normalization in ordinary space according to Eqs. (3.17) in 
Ref.0, $ffi(x, 0) = 1 and $^(0, &)rj= 1, do not carry over to the supersymmetric case. This is due the fact that the 
volume of some supergroups is zeroB resulting in the vanishing of $> kl2k2 (0, s ) f° r certain values of k\ and k 2 . This 



collides with the normalization of the plane waves (2.7) to unity at the origin. The reason of this contradiction is a 
well known phenomenon in superanalysis. In going from Cartesian to angle eigenvalue coordinates, one has to add 
additional terms to the measure to preserve the symmetries of the original integral. These are called Efetov-Wegner- 
Parisi-Souijlas terms in physical literature. A full-fledged mathematical theory of these boundary terms was given by 
RothsteirM 

To solve this normalization problem, we use the following strategy. First, we evaluate the supermatrix Bessel 
functions without taking care of the normalization. We just multiply the integrals with a normalization constant 
Gfc 1 2fe 2 - Having done the integrals, we determine the normalization by comparing the asymptotics of the supermatrix 
Bessel function for large arguments with the Gaussian integral. 



III. SUPERSYMMETRIC RECURSION FORMULA 



We extend the recursion formula in ordinary spaceB to superspac e. Aft er st ating the result in Sec. [II A, we present 



the derivation and the calculation of the invariant measure in Sees. IIIB and [II C, respectively. 



A. Statement of the Result 



Let$fe 1 2fc 2 (s, r) be defined through the group integral in Eq. (2.1). It has two diagonal matrices defined as in 
Eq. (2.2) as arguments. It can be calculated iteratively by the recursion formula 



$fe 1 2fe 2 (s,f) =G kl 2k 2 / dn(s',s) exp(i(trgs-trgs')rii) $ {kl -i) 2k2 {s' ,r) , 



(3.1) 



where <J>^ (s', r) is the group integral (2.1) over UOSv ((fci — l)/2fca)) and G kl2k2 is a normalization constant, see 



- (k x -\)2k 2 



the previous section and Sec. VI. As in the ordinary caseta, the coordinates s' are radial Gelfand-Tzetlin coordinates. 
Again, they are different from the angular Gelfand-Tzetlin coordinates which will be discussed elsewheret^l. We also 
introduced the diagonal matrix 



r = diag(r 2 i, . . . ,r klll ir 2 ) = diag (r i; ir 2 ) 
such that r = diag (ru,f) and the diagonal matrix 

s' = diag(s' n , . . ■,s[ kl _ 1)1 ,is' 2 ) = diag (si, is' 2 ) . 

The invariant measure reads 

dn(s',s) = 2 k z+ 1 ti B (s' 1 ,s 1 )fi F {s 2 ,s 2 )»BF{s',s)d[Z'}d[s' 1 ] 



(3.2) 
(3.3) 



M-b( s 1j s l) 

Vf(s' 2 ,s 2 ) 
Hbf(s', s) 



v4 



A1K) 



Up 2 =iU q 2 =i( is P2 - is' q2 y 

rfcl-l 



(3.4) 
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Here, we have introduced the differentials 



die] = n « . d ^ = n ■ ( 3 - 5 ) 
p— i p— i 

The domain of integration for the bosonic variables is compact and given by 

s P i < Spi - s (p+i)i > p = l,...,(ki -1) . (3.6) 
The fermionic eigenvalues is' p2 are related to Grassmann variables £' and through 

ICpP = *«p2 _ * S P2 ■ (3.7) 

The Jacobian or Berezinian consists of three parts. One of them, {1b(s[, s±), depends only on bosonic eigenvalues and 
one, iu,f(s'i, Si), only on fermionic eigenvalues, i. e. only on Grassmann variables. The third part mixes commuting 
and anticommuting integration variables. To underline once more the difference between radial and angular Gelfand- 



Tzetlin coordinates which is also present in superspace, we mention that the radial measure (3.4) is quite different 
from the angular oneE^I. 

As in ordinary space, the recursion formula is an exact map of the group integration onto an iteration exclusively 
in the radial space. Having done the iteration on the first k\ levels, we have treated all Grassmann variables. Thus, 
in the integrand, we are left with the matrix Bessel function <£ > jf 4 - ) (— i2s^ 1 1 \t%) for USp(2k2) in ordinary space0, 



. fci-i 



J n=l 

exp (tttag^-^-trg^W) ex P (^iTVi) Sff^s^-Va) . (3.8) 

We have set s — and s' = It is worthwhile to notice that the radial Gelfand-Tzetlin coordinates have a 

highly appreciated and valuable property: The Grassmann variables only appear as moduli squared in the integrand. 
Thus, the number of integrals over anticommuting variables is only half the number of the independent Grassmann 
variables. Moreover, advantageously, the exponential is a simple function of the integration variables. Thus, we may 
conclude that the radial Gelfand-Tzetlin coordinates are the natural coordinates of the matrix and the supermatrix 
Bessel functions, because their intrinsic features are reflected. 



B. Derivation 

All crucial steps needed for thcj-derivation of the supersymmetric recursion formula ([O]) carry over from the 
ordinary recursion formula in Ref.£3. We order the columns of the matrix u E UOSp(k\/2k2) in the form u = 
[u\ u% ■ ■ ■ Uk-t Ufci+i • • • Uk 1 +k 2 }- We also introduce a rectangular matrix b — [u2 ■ ■ ■ Uk t u^+i ■ ■ ■ u^+k-A such that 
u = [m b]. Analogously to the ordinary case, we have 

Wb = l(fe 1 _i)2fc a 

brf = u p u\ + u p v) p = lfc a 2ft 2 - u\u\ . (3.9) 

p=2 p=fei+l 

We define the square matrix a = b'sb and rewrite the trace in the exponent as 

trgu'sur — trgor + o\\r\\ , (3.10) 



with (Tn = u\sui. Similarly to the ordinary case, the first term of the right hand side of Eq. ( 3.10| ) depends on the 



last ki — 1 + &2 columns u p collected in b and the second term depends only on u\. Thus, it is useful to decompose 
the invariant measure, 

d/j(u) = d/J,(b) dfj,(u{) , (3-11) 



and to write Eq. (2.1) in the form 
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®ki2k 2 (s,r) = dfi(ui) exp(ia n r n ) / dp,(b) exp(itrgcrf) 



(3.12) 



Since the coordinates b are locally orthogonal to u\, the measure d[i(b) also depends on u\. 

We now generalize the radial Gelfand-Tzetlin coordinates introduced intij for the ordinary spaces to the superspace. 
Naturally the projector reads (lfe 1 2fe 2 — uiu\) and we have the defining equation 



(lfci2fe 2 - uiu{) s (l kl 2k 2 - uiu\) e' p 



s' e! 



1, . . . ,ki - 1, fci + 1, 



(3.13) 



for the (k\ — 1 + k 2 ) radial Gelfand-Tzetlin coordinates s' p and the corresponding vectors e' p as eigenvalues and 
eigenvectors of the matrix (lk 1 2k 2 ~ u i u \) s (lfe 1 2fe 2 — uiu\) which has the generalized rank k\ — 1 + k 2 . Due to 



u[e' p = 0, we find 



(lfci2fc 2 - uiu\) se' p = s p e p , p = 1, . . . , h - 1, k\ + 1, . 
As in RefR the eigenvalues s' are calculated from the characteristic function 

z(s' p ) = detg Ulfc l2 fc 2 - u-lu{)s - s' p ^ 
= -s p detg (s - s p ) u[ ui 



,k 2 



(3.14) 



(3.15) 



which has to be discussed in the limits 



for p = 1, . . . , k\ — 1 

oo for p = k\ + 1, . . . , k\ + ^2 



(3.16) 



Thus, together with the normali zatio n u\u\ — 1, these are k\ + ki equations for the elements of u\. 

The two parts of the integral ( 3.12 ) have to be expressed in terms of the radial Gelfand-Tzetlin coordinates s' p . In 
a calculation fully analogous to the ordinary case, we find 



a n = trgs - trgs' . 

The eigenvalues t p , p = 1, . . . , k± — 1, k± + 1, . . . , k± + &2 of a obtain from the characteristic function 
w(t p ) = detg (<t - t p ) = -— detg ((lfc l2 fc 2 - uxu\)s - t p ^j . 



(3.17) 



(3.18) 



Comparison with Eq. ( |3.15 ) shows that the characteristic functions w(t p ) and z(s p ) are, apart from the non-zero 
factor —t p , identical. This implies t p = s' p , p = 1, . . . , ki — 1, k± + 1, . . . , ki + &2- Thus, by introducing the square 
matrix u which diagonalizes a, we may write 



a = sb = u's'u . 

By construction, u must be in the group UOSp{k\ — I/2 A/2), because a and a share the same symmetries. 
These intermediate results allow us to transform Eq. (3.12) into 

®k 1 2k 2 ( s , r ) = / d^(s',s) exp(z(trgs - trgs')rn) / dfi(b) exp(itrguh'ur) 



(3.19) 



(3.20) 

where c£/i(s' , s) is, apart from phase angles, the invariant measure d/i(ui), expressed in the radial Gelfand-Tzetlin 
coordinates s'. To do the integration over b, we view, for the moment, the vector u\ as fixed and observe that the 
measure dfi(b) is the invariant measure of the group UOSp{k\ — 1/2/02) under the constraint that b is locally orthogonal 
to 1*1. The matrix u G UOSp{k\ — \j1k-i) is const ructed from b under the same constraint. Thus, since b and u cover 
the same manifold, the integral over b in Eq. ( |3.20 ) must yield the supermatrix Bessel function <f \kx-i)2k 2 ( s \ t) an d we 
arrive at the supersymmetric recursion formula (3.1). In the last step, we used a line of arguing slightly different from 
the derivation in ordinary space. In this way we avoided a discussion related to the ill-defined supergroup volume. 
The invariance of the measure is the crucial property we need for the proof and this holds both in superspace and in 
ordinary space. 
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C. Invariant Measure 



In order to evaluate the invariant measure, we have to solve the system of equations ( 3.15 ) for Ivp 1 ^ 2 



p = l,...,ki and \a p ^\ 2 = |u( fel+2 p)i | 2 + |u(fe 1 +2 P -i)i | 2 , P = 1, ■ ■ ■ ,k 2 in terms of the bosonic eigenvalues s p = s' pl , 
p = 1, . . . , k% — 1 and the fermionic eigenvalues s' ki+2p = s 'k. 1 +2p~i ~ * s j> 2 ' P = 1, ■ ■ ■ , fe- 



\Upl\ 



p=i 

fci |,,(1)|2 

^— I e_, — e' 
9=1 



Ei4 1} i 2 5 

P =i 



S gl S pl g=1 iS ?2 S p i 



p = 1, . . ., k\ - 1 



(3.21) 
(3.22) 

(3.23) 

In App. [A|, we sketch the solution of this system for small dimensions. Ins pired by t hese solutions one can conjecture 
the general solutions and verify them by plugging them directly into Eq. ( |3.21 ) to ( |3.23 ) , one finds 



„(!)|2 



>2 



k 2 

V- 



P 2 . 



00, p=l,. 



,(1)| 2 = 



p=l,...,k 1 , 



|a«| 2 = 2 K 2 - z^)^ 1 ^^ " < l)n *=^>* 2 " 



n g ll,g^p(* S p2 " «S^ 2 ) 2 ]lgll( is p2 - Sgl) 

p=l,...,fe 2 . (3.24) 

These expressions are reminiscent of the ones derived in Ref.0 for unitary matrices. However, importantly, all products 
in (3.24) involving fermionic eigenvalues are squared. This reflects the degeneracy of s in the fermion-fcrmion block. 
We have introduced a new anticommuting variables with |^| 2 = is' p2 — is p2 according to definition (3.7). 

From this point op-,, the invariant measure can be calculated in th e sa me way as for the angular Gelfand-Tzctlin 
coordinates, see Ref.ll2l for details. The result is summarized in Eqs. (3.4). 



IV. THE FUNCTION $ 22 (s,r) 



We use the recursion formula (3.1) to calculate the supermatrix Bessel function for UOSp(2/2). To avoid the 
imaginary unit in the exponent, we study $22 ( — is, r). The recursion formula reads 

(4.1) 



(4.2) 



$ 22 (-is,r) = G22 J dn(s',s) exp((trgs - trgs')m) $ 12 (-is',r) 

The function $>i 2 (— is',r) is easily found to be 

$i 2 (-is',r) = Gi2N--2(r 2 i - in 2 )(s' n - is' 12 ) J exp(2ri 2 si 2 ) . 

The measure of the coset UOSp{2/2)/ UOSp (1/2) is according to formula ( j3~4| ) given by 

{isxi - s' n ) Ul=i( is 'u - s m) 



d/j,(s , s) 



nLiKi - ( is i2 - s 'uf 



-ds' n d^d4[ 



We do the Grassmann integration and find 



$ 22 (-is,r) = G22 exp (r u (s u + s 2 i) - 2isi 2 ir 12 ) 

/S21 2 / 2 

Mb(s'js) 1^(^12 - s q i) I 4j^[(iri2 - r 3l ) 
0=1 ^ 7=1 



(4.3) 



G 



-2(ir 12 - rn) h 2Mn(s f 1 ,sx) J exp (s'nO^i - r n )) ds' n , 

q— 1 * ' 

(4.4) 



where we have introduced the operator 
M mj (s' 1 ,s i ) 



(is m2 — s'ji) 

I ki 1 fei 

ly * - y 



(4.5) 



V 



For later purposes, we introduced general indices m and j. Obviously, the Grassmann integration yielded eigenvalues 
in the denominator. This is somewhat surprising because of the following observation: we can always parametrize 
the group element u G UOSp(2/2) in a non-canonical coset parametrization in the spirit of an Euler parametrization 
in ordinary space. Inserting this parametrization into the defining equation of the supermatrix Bessel function (2.1) 
one can expand the trace in all Grassmann variables. The expansion coefficients are polynomials in the commuting 
integration variables and - more important - in the matrix elements of s and r. The invariant measure can be expanded 
in the Grassmann variables as well. It does not depend on r and s. Although this procedure becomes rapidly out of 
hand even for small groups, it is clear that the outcome of this expansion will be polynomial in the eigenvalues of s 
and r. In other words: eigenvalues can only appear in the denominator by an integration over commuting variables 
and never by a Grassmann integration. Therefore, before performing any integral over commuting variables, there 
must exist a form of < I > 22( — is, r), which is polynomial in the eigenvalues of s and r. 

To remove the denominators and to obtain such a polynomial expression, we use the following result: Let f(s'i) be 
an analytic, symmetric function 'm s' u ,i — \, . . .k\. Furthermore, define the operator 

L m ( S ) = X> . (4.6) 

.. \ lSm2 - Sjl dSjl 
3=1 J J 

Then the action of the operator on the integral over the bosonic part of the measure is given by 

L m (s) ... f n B (s', s) /K) d[ S [] = 

- f 21 ... f ^ fi B (s',s) AWs' llSl ) /(si) d[s[] . (4.7) 

This formula is derived in App. pi 

We now set /(s^) = exp(— s' 11 (r2i — rn)) and insert Eq. (4.7) into Eq. (4.4), we arrive at 

I 2 

$ 22 (-is,r) = G 2 2exp(-2isi 2 iri2) 4j^(ir i2 - rji)(is 12 - Sj±)- 

2~Yjisi2 - s q i) (iri 2 -r 2 i - rn - J^j^J ^ H^i, r i) , (4-8) 

where $2 1 '( s i: r i) is the matrix Bessel function of the orthogonal group 0(2) in ordinary space as defined in Ref.El. 
Although this can .already be taken as the result, we underline the symmetry between s and r by using the explicit 
form (3.21) of Ref IB for $ 2 (si,ri) 



$ 22 {-is,r) = G 22 cxp [trgrs - - 



2, 

2 2 2 



4jJ(w- 2 i - r lj )(is 2 i - sij) - y^X is 2i - siq) y (ir 2 i - r lp ) 

j=l q=l p=l 

- z -£j2nI (z/2) , (4.9) 
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where we have introduced z = (sn — Si2)(?"ii — m) and the modified Bessel function 1$ as defined in Ref.lil. 

The result fl4.7| ) was crucial in the derivation of $22 ( — is, r). By means of this formula, the denominator problem 
was overcome in one step. Because of its importance, we want to gain more insight into this problem: In App. ^|, we 
rederive $22 ( — is, r) in two other ways. It is clear that the methods of App. |c|can not be used for higher dimensions k\ 
and 2^2, but it will help to understand the mechanisms needed when working with radial Gelfand-Tzetlin coordinates. 



V. THE SERIES OF FUNCTIONS $ kl4 (s, r) 



We c alculate iteratively the four supermatrix Bessel functions <I>fc l4 (,s,r) for k\ = 1,2,3,4. We do this in Sees. V A 
to VD, respectively. 



A. First Level k\ = 1 



According to the recursion formulae (p.l|) and (3.8), the starting point is the matrix Bessel function for the unitary 



symplectic group < & 2 4 ' ) (s2, r 2 ), which was already calculated in Ref.Eii Up to a normalization, we have 

^ )(l2S2 ' r2) = |J 2 (a!(»* 2 ) A! (*(«*)) " A3(M 2 )A3( W (tr 2 ))) 

cxp(— tr is2w(ir 2 )) , (5.1) 



Since the subgroup 0( 1) o f UOSp(l/A) is trivial, no commuting integral has to be performed to derive $ 14 (— is, r). 
Inserting the measure (3.4) into the recursion formula and performing the Grassmann integrations yields straightfor- 
wardly 

3>i4(— is, r) — Gu exp (trg rs 



A 2 2 (ir 2 )A 2 2 (is 2 ) A 3 2 (ir 2 )A 3 2 (is 2 ) 
2(zs 2 i - sn)(ir 21 - rn) - 1 ) I 2(is 22 - sn)(ir 22 - r u ) - 1 



p- exp (trg (rs)) ,. , . „. 

A 3 2 {ir 2 )A 3 2 (is 2 ) 



The exchange term (ir 12 < — > ir 22 ) accounts for the permutation group S 2 in Eq. ( |5.l|) . Anticipating that the structure 
of $u(—is,r) will, remarkably, survive on a ll level s up to $ 44 (— is,r), we state that $ 14 (— is,r) essentially consists 
of two parts. A comparison with Eqs. ( |4.2| ), and (5.1) shows, that the first part of $i4(— is, r) is a product of an 
exponential with three other terms. The first one, 



1 



Ai{ir 2 )A 2 2 (is 2 ) A^(zr 2 )A^(z S 2) 



(5.3) 



stems from the integral over the USp(4) subgroup. The other two terms can be identified with the supermatrix Bessel 
functions 

$i 2 (-is,r) with s = diag (s n , isi 2 , is 12 ), r = diag (r n , ir\ 2 , ir 12 ) (5.4) 

and 

$i 2 (-is,r) with s = diag (su, is 22 , is 22 ), r = diag (rn, ir 22 , ir 22 ) . (5.5) 

The second part can be considered as a correction term, which destroys the product structure of <5>i 4 (— is, r). We may 
identify the different parts of the product with the integrations over the corresponding subsets of the group. Thus, 
$2 4 ^( — isi, t 2 ) arises from the integration over the USp(A) subgroup, the 0(1) integration yields unity and the other 
two factors come from the integration over the coset U0Sp(l/A)/(USp(4) <E> 0(1)). 
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B. Second Level ki — 2 



We now have to do one integration over a commuting variable. After the Grassmann integration, we are left with 
a considerable amount of terms. To arrange them in a convenient way we, introduce the following notation for the 
product of two operators Di(s)D 2 (s) acting on a function f(s), we define 



[D?{a) D 2 (s)} f(s) = D^s) D a (a) f(s) - (D^s) D 2 (s)) f(s) 



(5.6) 



This means, an operator with an arrow only acts on the terms outside the squared bracket. With this notation we 
can write 



$ 24 (-is,r) = G 24 exp(tr (r 2 s 2 ) + r 2 i(sn + s 2 i)) 

2 2 

duBis'n si) 



1 



i=ij=i 

i 



A 2 2 {ir 2 )A 2 2 (is 2 ) A 3 2 (ir 2 )A 3 2 {is 2 ) 

2 

r 2 i - «ri2 



\ i=l fe=l 
/ 2 2 

4jJ(ir 22 -r a )+2^ 

V i=l fc=l 



IS\ 2 - Sfcl 



r 2 \ — ir 22 
is 22 - Ski 



+ 2M^(s[,s 1 ) 



2M 2 i(si,si) 



A|(ir 2 )A|(i S2 ) A|(ir 2 )Af(£ S2 ) 



(- ^— rAf 21 (ai,si)- ^-Afn(ai,si) ) 

\IS12-S n IS 22 - S n J 



Al(ir 2 )A\{is 2 ) 
2 

'Al(*r 2 )A*(w a ) 



2tr r\ — tr ir 2 + 



i=l 



is 22 - Sn 



(o ^ 1 

ztrri — trzr 2 + > 



Afn(«i,Bi) 

M 2 i(s' 1 ,si) 



2 2 



A|(ir 2 )A|(is 2 ) ^ 1 = 1 Sfcl - is i2 



exp(s' u (ru - r 2i )) 



W" 22 ) 



(5.7) 



As in Sec. IV, a denominator problem occurs. It becomes obvious in the product M{^(s' 1; si)M 2 i(s' 1 , si). Thus, we 
expect an identity similar to formula (4.7). This identity should map a product of operators Li(s)L 2 (s) acting on the 
integral onto a product of operators Mn(s' 1 , si)M 2 i(s' 1 , si) acting under the integral. Neither the outer operators, 
L m (s), nor the inner ones, M m j(s), commute. Hence, the desired identity must be a non-trivial one. It is given by 
the following result. 



We have the same conditions as in formula (4/7), furthermore we define 



n— 1 q— 1 



l 



iSm% - S nl )[iSl2 - S q l) ds n ids q i 



(5.8) 



Then the following formula holds 



fiB(s',s)d[s' 1 ]f(s' 1 ) = 

s (fel-l)l 

fcl-1 fcl-1 

Vb(s',s) E M mj( s 'n s l) M lk(s'l,Sl)- 

3=1 k=l 
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fcl-1 



■E 



3=1 \ Ji 

fei— 1 



1 



— M m j(si,si) 



(is m2 ~ s' kl )(is m2 - s' ol ){isi2 - s' kl ){isi2 ~ s' ol ) 



M)d[ s [] . 



(5.9) 

The deri vation is along the same lines as the one for formula (|4.7|) , it also involves formula fl4.7|) . With the identi- 
ties (5.9) and ( |4.7D the denominator problem is again solved in one step. After some further manipulations we arrive 
at 



<& 24 (-is,r) = 27rG 24 exp (trgrs 



1 



A2(ir 2 )A2(is 2 ) A3(ir 2 )A3(i S2 ) 

2 

4 JJ(ra - iru)(sn - is 12 ) 

i=l 

- y^( s ji _ ^12) (ni - ^12) - z%— 

5=1 

2 

4 JJ(r»i - ir 22 )(sii - is 22 ) 



~ X^' 1 _ * s 2 2 )(^l - «>22) - Z — 



i=l 
3=1 



J («/2) 



-27rG 24 exp \trg rs - -J 

2 2 2 

A^A^) g Q( s « " " W2) 



-27rG 2 4 exp ( trg rs 



k=l ' 



) (trgrs - |J 
1 







A3(zr 2 )A3(i S2 ) 



((trg S )(trgr)-l)z— 7 (z/2) 
+ (iri2 < — ► «r 22 ) . 



(5.10) 



As in Sec. we used the composite variable z — (sn — si 2 )(rn — ri 2 ). A comparison with Eqs. (4.8) and (5.2) 
shows the similarity in the structures of < I ) 2 4(— is,r) and $i 4 (— is, r). The former also decomposes into two parts. 
The first part is a product, whose factors can be assigned to the integrations over the different submanifolds of the 
group in the same way as in the case of $14 (—is, r). The other one can be interpreted as a correction term due to the 
non-commutativity of the operators L m in formula (|5.9|). 



C. Third Level ki = 3 



This structure of <&kt4,(— is,r) emerging in the previous calculations is likely to be also present for arbitrary k\. 
However, for k\ > 2, we have so far not been able to treat the general case. Fortunately in important physics 
applications, one matrix argument of the supermatrix Bessel function has an additional twofold degeneracy in the 
boson-boson block. In this case, it is possible to carry on the recursion up to $44(— is, r) by extending the techniques 
developed for k\ = 1 and fc 2 = 2. Thus, from now on, we restrict ourselves to this case. 

At first sight, one might hope to achieve some simplification by applying the projection procedure onto the degen- 
erate matrix, because this results in a considerable simplification of the invariant measure. However, it turned out 
that the integrations are ea sier if one does the recursion with the non-degenerate coordinates. Hence, we use the 
measure as it stands in Eq. ( |3.4|) . We consider <£>34(— is, r) in the case that 

ri = diag (7-11, r 2 i, 7-21) . (5-H) 
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Having performed the Grassmann integral, one can arrange the terms in a way similar to Eq. (5.7). The complete 
expression and further details are given in App. We then can use formula (|5.9| ) and find after some further algebra 

$ 34 (— is,r) =4G 34 exp (trr 2 s 2 )(r 2 i - iri 2 )(r 21 - ir 22 ) 
1 1 



A 2 2 (ir 2 )A 2 2 (is 2 ) A 3 2 (ir 2 )A 3 2 (is 2 ) 

2 3 

4 Y\ ( r ki - iriz) W ( s kl - is iz) 

k=l k=l 

3 3 , 

+ 2 ^ Y\_ ( S J! ~ iSl2 ) ( ri1 + r21 ~ iri2 
2 3 

4 Y[ ( r ki ~ ir 22 ) Y\_ ( s ki - i« 22 ) 

k=l k=l 
3 3 

+2^ JJ(sji - is 22 j | /-,] - / ;,( - //22 

fc=l jjtk 

4 



X! II^' 1 ~ * S 12)( S J1 - * S 22) 



A 3 (ir 2 )A 3 (is 2 ) ^11 
r ii + r 2i + r n r 2i - (rn + r 2 i){ir 12 + ir 22 )+ 



ir\ 2 ir 22 - (rn + r 21 - ir 12 - ir 22 ) 



d 
dsn 



3 3 3 



A 3 {ir 2 )A\{is 2 ) 



Yl n^J 1 ~ is i2)(sn - is 22 ) 



i,k i=i i=i 



5 



9 



9sa Ssfei 



^(-iai.ri) 



+ (iri2 < — * ir 22 ) , 



(5.12) 



where ^^\si,r±) is the matrix Bessel function of the orthogonal group 0(3). We notice, that the structure of 
is, r) and $ 24 (— is,r) reappears in <& 3 4(— is, r). 



D. Fourth Level fei = 4 

In the calculation of <I>44(— is, r), we again consider the case that the matrix r is degenerate, 

ri = diag(rii,rn,r 2 i,r 2 i) . (5.13) 

The main problem is t o fin d a convenient representation for the matrix Bessel function ^g^-is^, ri) appearing on 



the third level in Eq. ( 5.12 ). It turns out that the representation derived in App. B of Ref.til is very well suited to 
our purpose. Due to the degeneracy in r\, the original threefold integral can be reduced to an integral over just one 
single variable 

= ex P (r 21 tr S ;) /•+" ^exp (.(m -^jf) 

vni ~ r 2i 7-oo nti V«ii ~ 



Here, we again neglected the normalization because we want to fix it afterwards as explained above. Similarly, 
'i 1 



$4 (— is\,ri) can be written as a double integral, 



(rn - r 2i )2 
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ah dt 2 \t 1 -t 2 \— -4 — —2 (5.15) 

lli=irin=lVSil-»*r> 

Singularities have to be taken care of appropriately. After inserting Eq. ( 5.10| ) into the recursion formula and per- 
forming the Grassmann integration, on e can arra nge the terms in a similar way as in the case of <1 > 34(— is, r). At 
this point, we notice that formulae (4.7) and (5.9) need to be supplemented by further identities. We state the most 
important one in the following. 



The same conditions as for formula (4.7) apply. Moreover, we define the operator 

- 4- 

2 



9=1 



1 



iSm2 - S q x ds 2 ql 



1 



( d 



d 



(iS m 2 — Sql)(Sql — Snl) \ds q l 8s n l 



(5.16) 



Then we have 



L m {s) 



f XB(s',s)d[s[}f(s' 1 ) 



d 



(5.17) 



HB (s',s) Y M mj (s'l,Sl) -r-y- 
s C=i-i)i \_j=i i 1 

Again, the proof is along the same lines as the proof of formula (4.7) and the proof of formula (5.9) in 

Thus, there is a family of rules to transform operators symmetric in sn L m (s),L m (s) acting onto an integral into 
an operator acting under the integral. We need one more such transformation rule which tells us how the product 



L m (s)^Li(s) 



transforms into operators acting under the integral. This formula and further details are given in 



App. ^. Collecting everything, we finally arrive at 

$44 (-is, r) = 4 G 44 exp (-tr (r 2 s 2 )) \\( r ii ~ ir2 ^ 



1 



1 



A 2 2 (ir 2 )Al(is 2 ) A 3 2 (ir 2 )A3(is 2 ) 
2 4 

8 IT( ri1 ~ iri2 ^ II _ isi2 ^ 
i=i j=i 
4 4 , 

+4>, If (sji - isu) ( rn + r 21 - ir 12 - 

i=l j^i V 
2 4 

8^Q(r 4 i - ir 22 ) \\{sji - is 22 ) 

i=l j=l 

4 4 , 

+4y^]~[(sji - is 22 ) ir n + r 21 - ir 22 - 

»=1 i^i 

16 44 

-A3(^)A3( <S2 )gn("ii-"«)(^- 

(r 



dsn 

d 
dsn 

- is 22 ) 



11 t r 2i + r n r 2i - (iria + ir 2 2)(rn + ^21) 

1 a x 

+ir X2 ir 22 + -trgr- — J 

8 4 4 4 

~ A3r \ a 4 / ■ TT^' 1 ~ * Sl2 ) II^ 1 ~ * S22 ) 

A|(ir 2 )A*(«s 2 ) ^JUL xx 



S 



^\-is u n) 

ir 2 2) ■ 



(5.18) 



We mention that in the derivation of this result we frequently used properties of the matrix Bessel functions $3 (si, ri) 
and $4 4 ^(si,ri) that only hold for the case that one matrix has an additional degeneracy. 
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VI. ASYMPTOTICS AND NORMALIZATION 



The asymptotic behavior of the supermatrix Bessel functions calculated in the previous section s is a usef ul check 
which als o allow s us to hx t he no rmalization constants. We find from the expressions in Eqs. (hi), ( |4.S| ) and in 



Eqs. (5.2), (5.1C), (5.12) and (5.1! 



hm $ Wfe2 (- !V ) = 2 W! G il2l 



Hi! IIm=l O'l - is m2) (m - ir m2 ) 

Al(is 2 )Al(ir 2 ) 



det [ex V {2si 2 r J 2)] l j=l k2 .lim, <f> ki >(-isi,n) 



(6.1) 



In the degenerate case, each degenerate eigenvalue contributes-.according to its multiplicity. The asymptotics of the 
matrix Bessel functions of the orthogonal group is given byLD'Ej 



lim* (1) f-i Sl /t n) -c^^^- 1 ^i/ 4 det [ exp ^" irml ^^ n - m=1 -- fcl 



where the constant can be found in Muirhead's boot 



Thus we find 



Pi(fci) _ r ( fc l/ 2 ) /s?/2-fei/4 

fci! 

^ kl2k2 (-is/t,r) = 2^^ t (( fe 1 -^) 2 +( fe 1 -^))/ 4 C^)G fcl2fe2 

det[exp(s n ir m i/t)] n ,m=i,...,fei det [exp(2s j2 rj 2 /t)] i):! - =1 fe2 
B klk2 (s)B klk2 (r) 



(6.2) 



(6.3) 



(6.4) 



for the asymptotic behavior. 

Orj_the other hand, the supermatrix Bessel function relates to the kernel of Dyson's Brownian Motion in super- 
spaceEH. Due to the normalization of the Gaussian integral, 



^-(( fcl - 2fe )2 +(fcl - 2fc2 ))/4 2 ^_ fa _ fcl/s 



<l[a]cxY> [ --( CT -p) ) = 1 



(6.5) 



the kernel 



/7T A -((fc 1 -2/c 2 ) 2 + (fc 1 -2fe 2 ))/4 



u£UOSp(k 1 /2k 2 ) 



(6.6) 



is also normalized. Since it is obviously connected with to the supermatrix Bessel function by 

£r ((fci " 2fc2)2+(fei " 2fe))/4 9^i-fe-^/2 



rW^,t) = (£) 



1 



2 « 2 - 



exp ( -- (trgs^ + trgr^) ) <$> kl2k2 (-is/t, r) . 



(6.7) 



we can fix the normalization by using the asymptotic behavior 

'7T\ -((fei-2fe 2 ) 2 + (fc 1 -2fe 2 ))/4 2 fe !- fc 2-fcl/2 



limr fclfe (s,r,t) = 



fci! fc 2 ! 

det [<5(sii - r J -i)] i [J . =1 fci det [5{s l2 - r j2 )] iJ=1 _ k2 
B kl k 2 (s)B kl k 2 (r) 



of the kernel. Comparing Eq. ( |6.4| ) with Eq. (6.8), we find 

23fc 2 (fc 2 -fci)+fc 2 /4-5fc 2 /2-fci/2 



G 



k t 2k 2 



7r ((fe 1 -2fe 2 ) 2 +2/c 2 -2 fc2 )/4 fc2!r(fci/2) 



(6.8) 



(6.9) 



We mention that this calculation also shows that the diffusion kernels of the one-point function and of the two-point 
function of Dyson's Brownian motiontll, i.e. the function T/ 2k \ k (s,r,t) which was denoted by T k (s,r,t) in Ref.til, 
indeed satisfy the proper initial condition. 
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VII. APPLICATIONS 



Although we focus in this contribution on the mathematical aspects, we now briefly comment on a particular kind 
of application. As the reader will realize, our results derived in the previous section are, in some sense, more general 
than what we need in those applications on which we focus here. We take this as an indication that explicit results 
for even more complex supermatrix Bessel functions can also be obtained. The results of the previous sections yield 
the kernels of the supersymmetric analogue of Dyson's Brownian Motion for the GOE and the GSE in the cases k = 1 
and k = 2. We do .pot present the physics background here. The reader interested in these applications is asked 
to consult Refs.QEdo for generalities and RefO, in particular Sec. 4.2, for the issue discussed here. In the present 
contribution, we use the same notations and conventions. We restrict ourselves to the transition towards the GOE 
and suppress the index c. The corresponding form ulae fo r th e tran sition towards the GSE are derived accordingly. 
We treat the one- and -two-point functions in Sees. VII A and VII B| , respectively. 

Forrester and NagaoQ derived expressions for generalized one-point functions of Dyson's Brownian motion model 
with Poissonian initial conditions. The used an expansions of the Green function in terms of Jack polynomials. Datta 
and Kunzu employed a supersymmetric technique to address the two-level correlation function of the Poisson GOE 
transition. They arrive at a finite number of ordinary and Grassmannian integrals which are still to be performed. In 
our approach, we also arrive at a representation of the correlation function in terms of a finite number of integrals. 
However, since we managed to integrate over almost all angular integrals in the previous sections, our result contains 
considerably less integrals, in particular, no Grassmannian ones. It has a clear structure due to the fact that, apart 
from two integrals, all others are eigenvalue integrals, i.e. live in the curved eigenvalue space of Dyson's Brownian 
motion. Moreover, since our formulae for the kernel are valid on all scales, our result is also exact for finite level 
number. 



A. Level density 



We use the result (|4.9|), derived in Sec. IV, for t he supermatrix Bessel function ^ 22 (— is, r). Using the replacement 
r — ► (x + J) and s — ► s/t and the relation ( |6.7| ), we obtain the diffusion kernel for the level density 




r 1 ( S ,x + J,<) = (2vr)- 1 /2^ 

(I 1 2 

exp ( — (sii -xi- J\) 2 - -(s 2 i -x\ - Ji) + -(isi2 -xi + Jif 

\ 

isiz - Sji) + yfosi2 - Sgl) . (7.1) 

9=1 / 

We take the derivative with respect to the source variable J\ and arrive at the level density 

((is 12 ~ Sll ) + {is 12 -s 21 ))B 21 ( S )z[ \ S )d[s} , (7.2) 

where the Berezinian is given by Eq. ( 2.10| ) for k\ = 2 and k 2 = 1. This result is exact for an arbitrary initial condition 
and for arbitrary N. In the case of a diagonal matrix H^ ' as the initial condition, we have 

J n=l E[j=l( S il + lE - HU) 1 ' 2 

and analogously for the GSE. 

In the limit t — > oo the stationary distribution of classical Gaussian random matrix theory is recovered. This can 
be seen by re-writing Eq. ([m]) for the rescaled energy x\ = Xi/t and the rescaled source variable J\ — J±/t, see 
alsoEil. In this limit the average over the initial condition yields unity and we arrive at an integral representation of 
the one-point correlation function of the pure GOE, 



Ri(xi) = j- 27r )3/2 S / exp (~( Sl1 ~ Xl ) 2 ~ ( S21 ~ Xl ) 2 + ( isi2 ~ Xl ) 2 ) 



14 



Sll - S 2 1 



1 



(lSi2 - Sn)(iSi2 - S 2 lJ \?Sl2 - S11 «Sl2 - «21 



(s n +ie) N / 2 (s 21 +ie) N / 2 



d[s] 



(7.4) 

where the symbol 3 denotes the imaginary part. Eq. ( |7.4| ) is equivalent to the classical expressions for the one-point 
functions as in Mehta's booktil. 

Finally, we state an integral expression for the one-point function for the case of a Poissonian initial conditions, see 
Eq. (5.1) of Ref.El. We have 



Z[°\s) = 



n 7 -=i(* s j2 - z) 



N 



(7.5) 



Inserting this initial condition into Eq. (7.2) yields the level density of a transition ensemble between Poisson regularity 
and GOE in terms of a fourfold integral. A further analysis will be published elsewhere. 



B. Two— point function 



The result ( |5lf ), derived in Sec. [v], for the supermatrix Bessel function <&44(— is, r) gives, with the replacement 
r — > (x + J) and s — ► s/t and according to Eq. (6.7), the diffusion kernel for the two-point function 



T k (s,x + J,t) = exp^-i(trgs 2 + trg(a; + J) 2 )^ <$> 44 (-2is/t, x + J) . 
The derivative with respect to the source terms leaves us with the two-point correlation function 

3 / ,n 28 ^44 ffsf w(0), ^ ^Ll(» s l2 - s k i)(is 2 2 ~ Sfci) 
R 2 (xi,x 2 ,t) = -jj^r I I B i2 {s)Z 2 '(a) 



(is 12 ~ is 22 ) 2 



1 



exp ( -- (trsf + 2x1 + 2 ^2 - 2 (« s i2 ~ xif ~ 2(is 2 2 - x 2 f) 



E 



1 



(isi2 - s k i)(is 22 - Sji) 
t 



X\ — t 



ds n 



x 2 — t 



d 



x\ — t 



ds k i 
d 



dsji 
d d 



(xi - x 2 ){isx 2 - is 22 ){isi 2 - s k i)(is 22 - 8ji) 

1 *3 

2 (xi - x 2 )(isi 2 - Ski){is 22 - Sji)(isi 2 - is 2 2) 2 \dsji ds k 
$i 1) (-2isi/t,a;) ) d[s] + (xi < — > x 2 ) 



x 2 — t 



ds k i 



(7.6) 



(7.7) 



The last line indicates that the inte gral with x± and X2 interchanged has to be added. This yields yet another 
simplification, since all terms in Eq. (7.7) antisymmetric under interchange of Xi and X2 drop out. We arrive at the 
expression 



R 2 (xi 



exp ( — - (trsj + 2x 1 + 2x 2 - 2{is\ 2 - Xi) - 2{is 22 - x 2 ) 



E 



Xi — t- 



d 



x 2 —t 







dski 



(isi 2 - Ski)(is 22 - sji) 

^ ) {-2is 1 /t,x)d[s] (7.8) 

The symbol denotes a certain linear combination of R 2 (x\, x 2l t) as explained in Ref.0. The normalization constant 
obtains from Eq. (|6.9[) and is given by G44 = 2(27r)~ 4 . This result is an exact expression for the two-point function 
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of Dyson's Brownian motion for every initial condition. Plugging in the initial condition of Eq. ( |7.5| ), we find an 
integral representati on of the two-point function for the transition towards the GOE. We notice that 2is\/t, x) 

is, according to Eq. ( 5.15 ), give»,as a double integral. 

As already discussed in Ref.tll, the kernels for the supersymmetric version of Dyson's Brownian motion are the 
same on all energy scales. Thus, the integral representation derived here for the two-level correlation function is, 
apart from the initial condition, the same on the so-called unfolded scale which-is relevant for physics applications. 
The initial condition on the unfolded scale is found along the lines given in Ref.EJ. 



VIII. SUMMARY AND CONCLUSION 

We extended the recursion formula of Ref.0 to superspace. Due to the group structures in superspace, we could 
restrict ourselves to the unitary orthosymplectic supergroup. As in the ordinary case, the recursion formula is an 
exact map of a group integration onto an iteration in the radial space. We used it to calculate explicit expression for 
certain supermatrix Bessel functions. i_. 

In ordinary space, we saw that the matrix Bessel functions are only special cases of the radial functionslij. We have 
not yet studied this further, but in our opinion a similar generalization is likely to also exist in superspace. 

It is a major advantage of the radial Gelfand-Tzetlin coordinates in superspace that the Grassmann variables 
appear only as moduli squared. Thus, the number of Grassmann integrals is a priori reduced by half. As we showed 
in detail, this is highly welcome feature for explicit calculations. As a remarkable consequence of this recursive way to 
proceed, the structure of the supermatrix Bessel functions is only very little influenced by the matrix dimension. We 
also saw that the basic structures of the supermatrix Bessel function for smaller matrix dimensions survive during the 
iteration to higher ones. The matrix Bessel functions in ordinary space show similar features. There, the structure 
of the matrix Bessel functions is much more influenced by the group parameter (3 than by the matrix dimension. 
However, as in ordinary space, it remains a challenge to find the structure of these functions for arbitrary matrix 
dimension. 

In interesting feature occured which sheds light on the general properties of the recursion. Total derivatives showed 
up in the integral over the commuting variables after having done the Grassmann integration. Since similar terms 
already occured in ordinary space, they are likely to be an intrinsic property of the recursion formula. Here, we 
succeeded in constructing a series of operator identities to remove them. This was a crucial step for the application 
of the recursion formula. A deeper understanding of these identities is highly desired. 

It should be emphasized that the total derivatives are no boundary terms in the sense of Rothstein. We showed 
in detail that such terms cannot occur because we always work in a compact space. Thus, according to a theorem 
due to BerezinH, the transformation of the invariant measure to our radial Gelfand-Tzetlin coordinates cannot yield 
Rothstein boundary terms. However, if further integration over the eigenvalues is required, such terms can emerge. 

As an application, we worked out some kernels for the supersymmetric analogue of Dyson's Brownian Motion. 

The radial Gelfand-Tzetlin coordinates are the natural coordinate system for the matrix Bessel functions in su- 
perspace. This parametrization represents the appropriate tool for the recursive integration of Grassmann variables. 
Once the particular features of this parametrization are better understood, they may allow for the evaluation of higher 
dimensional group integrals. 
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APPENDIX A: RADIAL GELFAND-TZETLIN COORDINATES FOR THE UNITARY 
ORTHOSYMPLECTIC GROUP UOSP(ki/2k 2 ) 

We wish to express the moduli squared of the elements of an orthogonal (k\/2k2) dimensional unit supervector 
in radial Gelfand-Tzetlin coordinates. To illustrate the mechanism, we start with the smallest non-trivial case, the 
group UOSp(2/4). We notice that there are at first sight minor, yet crucial, differences to the calculation in Rcf.Ej 
where we also started with the smallest non-trivial case. The set of solutions of the Gelfand-Tzetlin equations ( |3.2j| ) 
involves one bosonic and two fermionic eigenvalues. The eigenvalue equations read 
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9=1 \ S <?1 — s l 



g=l («S,2 



.(1)1 



,(1) 



E 



,(1)|2 



„(1)|2 



9=1 ysgi 



,(!) 



«S92 



-,(1) 



(A.l) 
(A.2) 
(A.3) 



where the last equation has to be solved in the limit z\ — > oo. The bosonic equation (|A.2|) has a unique solution 



s' n . Taking s' n as new parameter, Eqs. (A.l) and (A.2) can be solved, 



2 _ s pl ^11 
SpX Sq\ 



!-E ! 



lS k 2 - ggl l a .(!)|2 



p=l,2 



We insert these relations in Eq. (A.3) and obtain 



_ • (i)/ / _ • (i) 

Z 2 — 1 S 2 \ ^ 1 1 ^ ^ 2 



2 , . (IV / 2 

9=1 («S g2 T V fc=l * 



Cfc 



■ (1) 

«Sfe2 - 



i4 13 i 2 



with z\ — > oo. Here, we have defined the commuting variables 



]lg=l(^fc2 - Sgl) 



fc = 1,2 



(A.4) 



(A.5) 



(A.6) 



It remains to determine the set of solutions of the fermionic eigenvalue equation ( |A.5| ). To this end, both sides are 
inverted 



o = \{^ q 2 - *4 1} ) 2 ( i - E ~ — —) 



(1)|2 



9=1 



k=1 iS k 2 ~ «S 2 



211 



„(1)|2 



. (l)l"fe 
fe= l «S fc2 - ««2 



(A.7) 



We can now take the square root on both sides 



= IJ(is g 2 - 

9=1 



fe=l is fc2 ~ 



(1)|2 



n 



Ck 



„(1)|2 



4 ±± ■ (1) I"* 

fe=l lSk2 ~ *4 



(A.8) 



The most general form of the fermionic eigenvalue is 

2 



is 



(i) 



a + E a *=l a fe I 2 + ai2 n 



„(1)|2 



(A.9) 



k=l 



k=l 



After inserting this ansatz in Eq. ( A.8 ), we obtain two sets of solutions for the coefficients aio,an2 and a^- with 
i = l,2,j = 1,2 



IS 12 = IS\2 + Ci 



ClC 2 



IS 22 = IS22 + C 2 



IS12 - IS 2 2 
C\C 2 

is 2 2 - isi2 



G 



(1)|2 



(1)|2 



,.(1) |2 



a 



(1)|2 



(A.10) 



Remarkably, we have a\2 — a,2\ = 0. This allows us to write the nilpotent part of is' k2 as the modulus squared of a 
new anticommuting coordinate. 



is' k2 = is k2 + \£' k 



(A.ll) 



We solve the equations (A.1C) for lo^'l 2 , insert the results in Eq. (A.4) and arrive at 
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,(«)|2 



»l 2 



(ggl ~ gjl) IIn=l( S Pl ~ ^ S "2) S 



2 K2 



,l( s pl - 

(is P 2 



^2 



s'n)(«s P 2 - «s g2 )^ 



(is P 2 



* s g2) 2 nr, = l( is p2 ~ Snl) 



(A.12) 



The struc ture of Eq. (A.12) indicates the form of the solutions for groups of higher order as they were stated in 
Eq. ( 3.24 ). They are checked by inserting them directly into the Gelfand Tzetlin equations ( [3.23| ). The algebra 
needed is, although tedious, straightforward and similar to the one here. 



APPENDIX B: DERIVATION OF FORMULA 

The technique we use is an extension of the one developed in App. D of RefM First, we rewrite the integral in 
terms of functions. The left hand side reads 

L m (s) / (j,b(s',s) /(si) d[s[] Y[@(s k i - s' n ) Y[ 00a - ■ (B.l) 

J k<l Kn 

The integration domain is now the real axis for all variables. The action of L m (s) on the integral yields: 

// fei fei — 1 

k<l Kn J 

fcl 



fe<; 2<n 



We decompose the first term in partial fractions and find 

k\ fei — 1 



/ \M^s) EE??* vr m M) 



fci-1 



An integration by parts yields 

fei-i 



/ MB ( S ', S ) e (-^(si, *o) n e (*w - ^1) n - d t s 'i] + 

j=l k<l Kn 

m E — — + E • / ^ 

\ f-f is m 2 - s a dsn w m2 - s' ds' n 



(B.2) 



e • — — ?r- n e ^ - s u) n - s «i) d t s i] • (B-3) 
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n - s 'n) n - s «i) d ^ ■ (B-4) 

k<l Kn 

The derivatives of the functions yield S distributions. Upon integration of the S distribution the two terms in the 
last integral cancel each other. Hence the last term vanishes identically. This completes the proof. 



APPENDIX C: ALTERNATIVE DERIVATIONS OF $ 22 (s, r) 

We present two different alternative derivations. We do this in some detail, because the calculations give helpful 
informations on the role played by the radial Gelfand-Tzetlin coordinates. 

First, we use a angular parametrization of the coset UOSp(2/2)/ UOSp(l/2) by writing the first column of u E 
UOSp(2/2) as 



»1 



a\ 2 cos$\ 
/ 



yl — |a| 2 sin$ 
V2 a 



(C.l) 



This is jt-canonical way to parametrize the supersphere 5" 1 ' 2 that is isomorphic to the coset UOSp(2/2)/ UOSp(l/2), 
see Rcf.Eil. It coincides with a special choice of the angular Gelfand-Tzetlin coordinates. The invariant measure-is in 
these coordinates simply dfi(ui) = da* dadd . Thus, one directly obtains the the volume V(S^ 2 ) = 0, see RefH In 
the parametrization of the measu re b y radial Gelfand-Tzetlin coordinates (4.3), one has to perform the Grassmann 
integration and to apply formula (|4.7|) to achieve this result. 



Although we use a diffe rent coordinate system, we still take advantage of the r ecur sion f ormula (3.1). To use it in 
the parametrization (C.l), one has to so lve t he Gelfand-Tzetlin equations ( 3.21 ) to ( |3.23 ) for the eigenvalues. The 
unique solution of the bosonic equation (3.22) is 



a 



The fermionic equation yields 



ll;=i0ii ~ ao ) | a |2 
isi2 - a 



a 



Sll + s 2 i Sn - S12 



cos2z9 



(C.2) 



is' 12 = is 12 + -\a\ 2 

IS12 - a 



(C.3) 



After inserting Eqs. ( |C.2| ) and ( |C.3[ ) and the measure dfi(ui) into the recursion formula (3.1), the Grassmann inte- 
gration can be performed. Remarkably, we arrive at the denominator-free expression 



2^ . 

$22(5, r ) — G22 I d-d exp ( trg rs — - + - cos 2§ 



2 1 2 

Y\_(rn ~ iri2)(sn - is 12 ) + - }Xsji - isi 2 )(ni - ir 12 ) 



irii - -(m +r 2i ) (sn - s 2 i)cos2i? 



; (ir 12 - r 2 i)(sn - s 2 i) sin 2 2^ 



(C.4) 



To make contact with Eq. ( ]4.9| ) one has to realize, that in Eq. ( p.4| ) an additional total derivativ e app ears in the 
integrand. This becomes obvious if one adds and subtracts z/4cos2?9 in the squared bracket of Eq. (C.4) 



$22(s,r) = G 2 



(z z 
trg rs — — + — cos 2t9 



Y\_(rn - ir 12 ){s l i - is 12 ) 
i=l 



19 



1 2 z 
-p / foji -*si2)(^i -iri 2 )- -cos2i? 



3=1 

W"i2 - r 2 i 



rai - 7-n Jq 



dx) 



d' 2 



o exp ( trg rs — - (1 — cos 2$) 
d(2tf) 2 V s 2 V > 



(C.5) 



While the first integral reproduces Eq. (4.9), the second one vanishes identically. In general, in performing Grassmann 
integrations, one has to take care of boundary contributionsaO. These contributions can appear whenever even 
coordinates are shifted by nilpotents and the function one integrates does not have compact supportcl. However, in 
our case the base space is always given by a n dimensional sphere, i. e. by a compact manifold wit hou t boundary. Thus 
in a properly chosen coordinate system, no boundary terms should appear. With regard to Eq. (C.5) this means: the 
fact that the last term vanishes, is a direct consequence of the compactness of the circle and of the analyticity of the 
function, that we integrate. However, in the radial Gelfand-Tzetlin coordinates, only the moduli squared of th e vector 
Mi are determined. Therefore, not the whole sphere, but only a (2 ra+1 ) th segment of it is covered by Eq. ( 3.24 ). In our 
case, not the circle but only a quarter of it is parametrized. This is allowed since the supermatrix Bessel functions 
depend only on the moduli squared |u i:L | 2 . Nevertheless, one has to ensure that the introduction of these artificial 
boundaries does not alter the result. To this end we use the following integration formula. 
Let sn < s' n < S21 be real and let £',£'* be anticommuting. Furthermore, define 



with two analytic functions /o(s'n)) /i( s 'n)- Then the integral 

Jsn 

transforms under a shift of s' n by nilpotents 



'ii 



in the following way, 



/ = 



S21 a I 

dy^*d^-^f(y(s' 11 )^,e)-[fo(s 2 i)g(s2i)-fo(su)9(su)} 



(C.6) 



(C.7) 



(C.8) 



(C.9) 



The proof is by direct calculation. The second term in Eq. ( |C.9] ) is often referred to as boundary term. It can be 
viewed as the integral of a total derivatbfe, i.e. an exact one-form, that has to be added to the integration measure 
for functions with non-compact supported. For functions of an arbitrary number of commuting and anticommuting 
arguments, a similar integral formula holds with additional boundary termsa. In going from the canonical coordinates 
($,a, a*) t o the radial ones (s'n, m principle boundary terms can arise, since t he b osonic Gelfand-Tzetlin 

eigenvalue (C.2) contains nilpotents. However, the crucial quantity is g(y) in formula flC.9| ) which, in our case, is 
given by 



isu - s'n 



(C.IO) 



Thus , g(s \ ; ) causes the boundary term to vanish at su and s 2 \- It is the product structure of the left hand side of 
Eq. ( [3.24 ) which always guarantees the vanishing of the boundary terms, when one goes from the Cartesian set of 
coordinates to the radial Gelfand-Tzetlin coordinates. 

Therefore, one may think of the denominators, arising in Eqs. (4.4), (4.5), as belonging to total derivatives of 
functions, which vanish at the boundaries. Keeping this in mind we derive Eq. (O) in yet another way. We expand 
the product 



W (is m 2 

9=1 



Sql) 



kl x 

71 = 



d" 



n 

9=1 



(Gil) 



and insert it into the integral 
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Sll 

. . . / MB(si) S l) IJ ( is rn2 - S„l) M,^^, 

-« JlVn-i)! n =l 



s 1 )K mj (s' 1 , si)/(s' 1 )d[s' 1 ] 



S(fc!-l)l 
Sfc, 1 



si)MM*'i] 



(C.12) 



We can remove the term proportional to {is r 



in the integrand by an integration by parts. Through the 



expansion (C.ll), the vanishing of the boundary terms is assured. We arrive at 
i^ TO j(si,Si) = - ^ — (is m2 - s^)"" 

n=2 

rigil ( !S «i2 - Sgl) ~ UgLl ( S jl - Sql) 



fcl -, ^„ fcl 

n=2 °V S ilJ g=1 

rfci / • _ \ rrfci 



S J1 



fcl 



fcl 



ly, 1 1 y 1 y 1 

^ 0=1 * Sm2 _ S 9 1 ^ Sji — S g i 9=1 Sj 1 — 



We notice that in the new operator K m j{s' 11 s\) all denominators of the type (is m 2 
fci = 2, we calculate 

Ku = - (isi2 + s'n - sn - S21) 



9s 



gi 



which can be inserted into Eq. (4.4) by using the definition ( C.12 ). Finally, the result 
substitution 



(C.13) 

-1 have disappeared. For 



(C.14) 
) is reproduced by the 



sn + s 2 i sn - S12 



cos2?9 



(C.15) 



see Eq. (C.2). In other words, we have seen that the result of this procedure is summarized in formula (4.7). 

Finally, some remarks are in order: First, from this discussion, one might conclud e tha t the radial Gelfand-Tzetlin 
coordinates are less adapted to the problem than the canonical parametrization (C.l), because, in the latter, no 
denominators appear. We stress that this is not true. Certainly, the denominators appear due to the shift of the 
bosonic variable by nilpotents in Eq. ( |C.2| ). However, the difficulty in deriving Eq. (4.S) is the identification of 
the different parts of the integrand after the Grassmann integration. Some of them belong to total derivatives and 
this problem exists in both parametrizations. Second, we emphasize, that the appearance of total derivatives in the 
integrand is not a peculiarity of supersymmetry. Already in Ref.Ej where the matrix Bessel functions in ordinary space 
were treated we had to solve a similar problem. The appearance of these total derivatives is an intrinsic property of 
the recursion formula. A geometrical interpretation of this phenomenon is highly desired. 



APPENDIX D: DETAILS FOR THE DERIVATION OF <f> 34 (~is,r) 

We always consider the case that one matrix has an additional degeneracy according to Eq. ( 5.11 ) and ( 5.13] ) . We 
introduce the notation 



Si 



(Si 



is 7 2) and i? 4 



mi - ir j2 ) 



(D.l) 



Due to the degeneracy, < &24(~ is,r) simplifies enormously as compared to the general result ( [5.10| ). We insert it into 
the recursion formula (3.1) and do the trivial int egra l over the 0(2) subgroup. After performing the Grassmann 
integrals we arrive at an expression similar to Eq. (5.7), 



®3i(—is, r) = 4 G 3 4exp (tr r 2 s 2 + m(sn + s 2 i)) 



/ d£tB(s'i,Si) JJ-Rii 

4=1 j = l 



21 



A 2 (ir 2 )A 2 (is 2 ) A3(ir 2 )A3(is 2 ) 



4 n - 2 e R ^ s ki + 2 e M y (*i . s i ) 



j=i 

2 



fe=l 
3 



2 



4 n r ^ - 2 e ^ ^ + 2 e m ^ ' si ) 

i=l fc=l i=l 

1 1 

A 2 (zr 2 )A3(z S2 ) + Ai(ir 2 )A4(i S2 ) 



E 7 



is 12 - s 



—M 2j {s / 1 ,s 1 ) - 



is 22 - s 



— Mij(si,si) 



n 



A 2 (zr 2 )A3(zs 2 ) A3(ir 2 )A3( lS2 ); 11 is k2 - s'^ 
2 



2 



A|(ir 2 )Af(i S2 ) 
2 

Ai(*r2)A*(w 2 ) 



V i=l / j=l 

/ 2 \ 2 

trgr + m - ^ ^i 1 E M M> Sl ) 



3 2 



Ai(<»>)A8«ft)Sn j W 

exp((s 21 + s' n )(r 21 -r u )) 
+ (iri2 < — ► ir 22 ) . 



(D.2) 



Formulae ( |5.9| ) and (4.7) are needed to remove the denominators, in a way similar as for $ 2 4(s,r). A single sum 
Ylj=i Mij(s\, s i) transforms according to formula (4.7). Moreover, we observe that parts of Eq. ( p.2| ) together 
with the product J2j=i ^-ij{ s 'ii s i) Sfe=i ^2k{s' 1 , s%) yield exactly the integrand of formula (5.S). Thus, it can be 
transformed accordingly. After rearranging terms, we arrive at the result (5.12). 



APPENDIX E: DETAILS FOR THE DERIVATION OF $ 44 ( -is,r) 



For the recursion, we need < &34(s, r) with degenerate r = diag (ru, r 2 i, r 2 i) according to Eq. (5.11). Using the 



integral representation (5.14) for $3 (— is'^ri) we find the helpful identity 



d d . (1) . „ 



1 



1 



2 ^ - ^1 



g 



d 



ds 



exp (-r 2 itra' 1 ) ^\-is[, n) 



(E.l) 



J 1 . 



We stress that this relation, which is crucial in the derivation, only holds, because of the degeneracy in the matrix r\. 
Employing Eq. (5.14) and another identity, 







E 7j7~ CXP (~ r2ltr s 'i) $ 3 H s 'i, ri) 

4 — 1 * 1 



(ru - r 2i ) exp (-r 2X tr si) $^ (-is[, n) 



(E.2) 



We insert $34(5',?) into the recursion formula (3.1). We then can arrange the terms emerging from the Grassmann 
integration in a way similar to the former cases. We obtain 

. 2 4 

$44(-«s,r) = 4 G 44 exp(trr 2 s 2 + rntrsi) / dfi B {s'i, si) R21 S jt . 

i=l j=l 



22 



A 2 (ir 2 )A 2 {is 2 ) A 3 2 (ir 2 )A^(is 2 ) 



8R 21 R 2 n - iR n R 21 Ski + 4 E [R21 

' 4 3/ 

%R 22 R\ 2 — ARi 2 R 22 S' fc2 1 + 4 I i?22 



9s' 



fe=i 



M 2i (si,si) 



»"21 - Til 



[ R 22 R\ 2 + R11R12 + R11R21 + 2^ 12 ^ 22 ) ^ii 1 J J 



A|(ir 2 )Af(i S2 ) £; 
9 



r 2 i - rn 



R11R21 + R11R12 + R12R22 + 2^ 12 ~*~ ^ 22 ^ E ^s^ 1 ^ J 



fc=l i,j 



exp (-rntrsi) ®jp (-is[, n) 



+ C(s,r) + (ir i2 < — ► ir 22 ) 



(E.3) 



Again, all operators with an arrow are understood to act only onto the term outside the squared bracket, i.e. onto 
exp (— rntr s[) (— is'i, Fi). In the function C(s,r), we summarized the terms that are expected to arise due to 
non-commutativity o f so me operators acting on the integral and some operators acting under the integral. The last 
two lines in formula (5.E) are examples of such terms 



. 2 4 

C(s, r) = 4 G 44 exp (tr r 2 s 2 + r u tr si) / dfj, B (s'i, s%) R 2l S 

J A 1 „- 1 



i=l j=l 



A|.(ir 2 )A|(i8 2 ) Ai(zr 2 )Al( ?S2 ) 



Y RnRu + (Ru + Ri2)(r 2 i - m) - (Ru + ^12)^-7 



^812 -sji «s 22 -s^i ^ 

1 



A^ir^A 2 .^) A3(zr 2 )A3(is 2 ) 



2 3/ \ 

n n Eni?i2 + + i? i2)( r '2i - ru) - (i?n + fli 2 )^r- 

fc=ii=i\ 



is fc2 - s 



1 



1 



A 2 {ir 2 )A 2 {is 2 ) A3(zr 2 )A3(i S2 ) 
^((r-2i-rn N 



ds 



(rax - r u ) - ) M^Ma* 



A'l 



23 



At(ir 2 )A 3 2 (is 2 ) 



da'. 



A|(ir 2 )A3(is 2 ) 



' jl / \fc=l 

9 .9 



J 2 



fc=l 



J ' 2 I ^1 



cxp (-mtrsi) $p(-isi,ri) 



(E.4) 



In order to evaluate Eqs. ( E.3 ) and (E.4) we need some more properties of the matrix Bessel function isi, ri). 

We investigate the action of Lk on $^(— isx,ri) using the integral representation ( 5.15 ). 
After a straightforward calculation involving an integration by parts we find 



L k exp (-r u tr s x ) $ 4 (-is lt r x ) 



± — 



i=l 



»Sfc2 - sn 



(r 2 i - r n ) 2 + (r 2 i - ru) 



9si 



exp (-rutrsi) $4 (-isi,ri) 



Now Eqs. flE.3|) and (E.4) can be enormously simplified by the observation that 

(r 2 i - r n )L k - L k ) exp(-r n trs 1 )$ 4 1) (— isi,ri) = 



which follows directly from Eq. (E.5). We find for Eq. (E.3) 

r 2 4 

®u(.-is,r) = 4 G 44 exp (trr 2 s 2 + rutrsi) / dn B (s[, si) 5, 



i=i fe=i 



A2(*r 2 )A2(* S2 ) A3(ir 2 )A3(is 2 ) 



R12R11 — 4-R 22 ^ 

4 3 

fe=i 3=1 

R12R11 I 8i?ni?2i — 4i? 2 i S^ 1 J 



fe=i 



8i? 22 i? 12 - 4i? 22 ^ S^ 1 + 4 ^ M 2j (si , Sl ) 



fe=i 



+ E i?11 ( r 2i - rn - — ) M^{s' 1 ,s 1 )M 2j ( S ' 1 ,s 1 ) 



J-' 



(r 21 -r u - —J Mi^^i.sijMjy^.ai) 



+ R2iR 22 M^ (a' x , si)M 2 i(s'i, si) 

A|fa 2 )Af(i S2 ) fliifii2 (^-JX 1 ) E( M «(«i.-o-^(«i.-o) 
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Al(ir 2 )Al(i S2 )En^' 5 fc/ 



exp (-rutr si) (-isi, n) 



fc=l i,j 

+ C(s,r) + (ir 12 < ► ir 22 ) 



The terms contained in Eq. (E.4) simplify, too. We arrive at 



(E.5) 



(E.6) 



(E.7) 
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r 2 4 

C(s,r) = 4 G 44 cxp (tr r 2 s 2 + rntrsi) / dlu, B (s[, s\) JJ Ru JJ Sji 

J „-_i 



«=i i=i 



+ 



A5(ir 2 )Al(w a ) A|(ir 2 )Af(Ma) 

3 f r 
R nRi2 + (^11 + Ri2)(r 2 i - r u ) - (i?n + #12)7^7 



i=i 



ii . 



16 



1G 



\ IS12 - Sjj 

1 1 



—M 2 j(s' 1 ,s 1 ) - — -j-Mijis'^si) 



IS22 - s 



A 2 2 (ir 2 )A 2 2 (is 2 ) Al{ir 2 )Al(is 2 ) 



! J i?lli?l2 + (^11 + fii 2 )(r 2 i - m) - (fin + flia)^7 



fc=l 

3 = 1 



<9~ 



9Sjl / ? ' Sfe2 _ S 'i 



exp (-mtrsi) $3 (-«'si,n) 



(E.8) 



To further evaluate the expressions, we can now invoke a symmetry argument between the eigenvalues ru and r 2 i, 
respectively. Since the product i?ni? 12 appears as a prefact or in front of t he integral (E.7), R 2 \R 22 must also appear 
as a prefactor in the final result. Thus, all terms in Eqs. (E.7) and (E.8) which do not contain R 2 \R 22 as a f acto r 
must yie ld zero. The remaining terms which are proportional to R 2 \R 22 can again be treated using formulae (4.7) 
and (5.£). However, we want to show explicitly that this line of arguing is correct and that the other terms indeed 
vanish. To this end, we need an additional identity to treat the operator product 



2 2 
I]^ M i7( s i> s i)^ M2fe ( s i> Sl ) 

3 = 1 ^J" 1 fc=l 



(E.9) 



The required identity is given by the following formula: The same conditions as for formula (4.7) apply, furthermore 
we define 



L Zi s ) L n{s) = 



1 



d 3 



-? (is m 2 - Sii)(is n2 - Sji) dsnds% 



(is m 2 - sa)(is n2 - Sji) 



d X 1_ 



dsn ~^ Sji - Ski \dsjx ds kl 



d d 



(E.10) 



Then we have 



L m( s ) L n(s) 



S21 /" s <=ii 



s-21 r s kll 



f x B (s',s)d[s' 1 }f(s' 1 ) 



iS n2 - lS m2 



(fcl-l)l 
fei-1 

E 

i=l 



(fci-i)i 
fci-1 fei-1 



3=1 i=i 



1 

— -^-T M mi[Sl,Sl) 



-y 

2 ^ (lSm2 ~ s' kl ){lS n2 - S' kl )(s' kl - S' n ) 2 8s kl 



4e- 



2 k^l ( iSm2 ~ S fcl)(* S ™ 2 ~ S il)( S fel ~ s ll) 2 ds kl 



f(s' 1 )fj. B (s',s)d[s' 1 ] 
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(E.ll) 

The proof is similar to the one of formula ( f4.7| ). We notice that the arrow in Eq. ( E.lCp is used slightly differently than 
previously. The operator L^(s) acts also on a part of L n (s). This is not consistent with the definition in Eq. fl5.6|) . 
However, since this is obvious where it occurs, we still use the same arrow. We can now translate the left hand side 
of Eq. (E.12 ) int o an e xpre ssion in terms of (— is, r). After some further manipulations involving the identities in 
Eqs. (E.6), (E.l) and (EJ2) we arrive at 



<^a{-is,r) = G44 exp (trr 2 ,s 2 + rntr s\) JJ-Rji \\ Ski 

i,j k—1 

\Al(ir 2 )A 2 2 (is 2 ) + A|(ir 2 )Al(i S2 )) 

^8i? 12 i? 22 - 4i? 22 J2 S k2 - 4£ 2 (s)^J \ 8R u R 2 i - AR 2l £ Srf - iL^sU 



A|( ? r 2 )A4( lS2 ) 
1G 



4 2 



A^r 2 )A^(is 2 ) 



j2 n R v s kj 



k=lj=l 

(ir 12 < — ► ir 22 ) 



cxp(-rntrs' 1 )$^ 1) (s' 1 ,ri) 



(E.12) 



After rearranging terms this yields the result ( 5.18 ) for $44 (— is, r). 
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